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! Abstract 

I We classify gamma-ray bursts (GRBs) according to their observed durations and 

physical properties of their spectra. We find that long/hard bursts (of duration Tqq > 2.5 s, 
> ■ 

' and typical photon energy Ep 0.8 MeV corresponding to BATSE's energy fluence 
hardness > 3) show the strongest deviation from the three-dimensional Euclidean 
brightness distribution. The majority of GRBs, i.e., short bursts (Tgo < 2.5 s) and long/soft 
bursts (with Tgo > 2.5 s, and H^2 < 3) show little, if any, deviations from the Euclidean 
distribution. These results contradict the prediction of extragalactic GRB models that 
the most distant bursts should be the most affected by cosmological energy redshift and 
time-dilation (long/soft GRBs). The strongly non- Euclidean GRB subclass has very hard 



' spectra of typical photon energy above ~ 1 MeV, i.e., outside the ideal energy range for 
I optimal detection by BATSE. Selection effects due to BATSE inefficient triggering and/or 
intrinsic burst properties may cause the strong deviation from the Euclidean distribution for 
the hardest class of long GRBs. Our results imply strong luminosity and spectral evolution 
of GRB sources placed at extragalactic distances. Whether there co-exist multiple GRB 
populations with different (Galactic or extragalactic) origins remains an open question. 
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1. Introduction 

The brightness distribution of gamma-ray bursts (GRBs) depends on both the intrinsic 

properties of the emission and the spatial distribution of sources. When analyzed globally, 

GRBs detected by BATSE show a clear deviation from the three-dimensional Euclidean 

brightness distribution (Meegan et al. 1992, 1995). However, recent analyses of special 

GRB subclasses showed that the deviation from the Euclidean distribution does not 

apply to all bursts. On the contrary, 'soft' bursts with typical long durations appear 

to have an Euclidean brightness distribution (Pizzichini 1995, Kouveliotou 1996, Belli 

1997, Pendleton et al. 1997). A commonly used spectral parameter is the hurst-averaged 

hardness ratio Hz2-, i.e., the ratio of overall fiuences in the BATSE energy channels no. 2 

and 3 corresponding to the 100-300 keV and 50-100 keV bands, respectively. The photon 

fluence-averaged ratio H^2-i ^^'^ energy fluence- averaged ratio i7|2 are commonly used. 

A first approach was to consider different GRB subclasses ranked according to their 

(e.g., Pizzichini 1995). The analysis of Kouvehotou et al. (1996) found that GRBs with 

photon-fluence-averaged > 4 and i?32 < 2 are distributed quite differently, with the soft 

bursts showing an Euclidean distribution. Belli (1997) used a different phenomenological 

criterion to distinguish GRB classes in the hardness/duration diagram {H^2 ^90) with 

T90 the duration of 90% energy deposition in the 50-300 kcV band) divided by the relation 
1/2 

H^2 — 229Q . The long-duration subclass with H^2 > 4.5 shows the most pronounced 
deviation from the Euclidean distribution. 

Pendleton et al. (1997, hereafter P97) used a very different approach reaching 
qualitatively similar conclusions. P97 select GRB individual pulses according to the 
existence or lack of a high energy component above 300 keV (BATSE's channel no. 4). Two 
classes of GRBs with without a high energy component (HE and NHE) were identified. 
Very often, complex GRBs with multiple pulses show both kinds of HE and NHE pulses. 
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The interesting conclusion of this analysis is that NHE pulses, both individual or as part of 
more complex GRBs, appear to have Euclidean distributions (P97). 

All of these analyses were carried out on a purely phenomenological level, with no direct 
correspondence between the selecting criterion and a physical model of emission. The main 
motivation of our work is to provide a physical model for the GRB observable quantities 
(T9o,-ff32) to better understand the GRB brightness distribution. We use the shock 
synchrotron model (SSM) of GRB emission that successfully explains the great majority of 
GRB burst-avcragcd spectra (Tavani 1996a,b, 1997; hereafter T96a,b, T97). The results 
of our GRB analysis are obviously independent of the SSM: only the interpretation of the 
adopted GRB selection depends on the theoretical model. We studied GRBs of the Fourth 
BATSE (4B) Catalog (Paciesas et al. 1997). Three effects should be considered: (1) the 
intrinsic spectral evolution within bursts (typically hard-to-soft) that tends to concentrate 
bursts diagonally from top left to bottom right of the H^2 ~ ^90 diagram; (2) the effect due 
to cosmological redshift, or any other cooling mechanism simultaneously dilating the time 
profile and softening the spectrum (T96a,b); (3) selection effects due to BATSE sensitivity. 
If the majority of GRB sources are at extragalactic distances, an overall softnening of 
dim and time-dilated bursts is expected as these bursts should be the most distant for no 
spectral cosmological evolution. The population of soft and long duration GRBs is expected 
be inhomogeneously distributed if BATSE's sampling radius exceeds the spatial distance of 
peak burst formation rate. As we show below, this is the opposite of what we find. 

2. Physical interpretation of Tqq and Hs2 

The double-peaked distribution of GRB durations is well known (Fishman & Meegan 
1995), and probably reflects different manifestations of simple (short) and complex (long) 
burst emission. Pulse analysis of complex bursts shows a structure of pulses with average 
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durations of ~ 200 — 400 ms (Norris et al. 1996). The subclass of short bursts with 
Tgo < 2.5 s can be modelled as relatively simple bursts whose substructure consists of a 
single (or a few) pulse(s) of average properties. Longer bursts (Tgo > 2.5 s) are either 
apparently dilated versions of simple short bursts, or complex superpositions of many pulses 
of average properties. The double-peaked Tgo distribution can then be related to the relative 
complexity of the pulse structure. No apparent correlation exists between overall burst 
durations and peak photon intensities. We define here 'short' and 'long' bursts divided by 
the critical duration Tgo = 2.5 s corresponding to the local minimum of the Tgo distribution. 

There is instead a definite correlation between Tgo and average GRB hardness, most 
hkely as a consequence of the hard-to-soft spectral evolution within bursts (e.g.. Ford et 
al., 1995). Fig. 1 shows the 4B hardness/duration diagram (H^2 vs. Tgo) with a clear 
indication that short bursts are generally harder than long ones. Spectral softnening for 
individual pulses and for long bursts has been reported in many cases (Ford et al 1995), and 
it constitutes a fundamental property of cooling particle distributions (e.g., Tavani, 1998). 

We apply the synchrotron shock model (T96a,b) for the interpretation of the hardness 
ratio H^2- The SSM is based on a hydrodynamic and radiative model of relativistic 
MHD outflows interacting with shocks that reconvert part of their energy in impulsive 
acceleration. The SSM can in principle be applied to relativistic mass/electromagnetic 
outflow models at extragalactic and Galactic distances. Average hardness ratios H^2 can be 
reliably converted into average peak energies, Ep, of the u F^, spectrum (T97). Fig. 3 of 
T97 shows the SSM calculated H^2 (and iffs) as a function of Ep for different assumptions 
about the index S of the high-energy component of the post-acceleration particle energy 
distribution function {N{'y) oc 7"''). The SSM calculated H^2^s are in good agreement 
with detailed broad-band average spectra of GRBs (T96a,b;T97). The agreement between 
theory and observations extends also to GRBs detected above 1 McV (T96a,b). a factor of 
~ 2 with the SSM calculated value an apparent depletion of soft photons We can therefore 
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physically interpret the average hardness ratio by converting its value to a typical photon 
energy Ep at which most of the burst energy is dissipated. We choose H^2 = 3 as the 
critical value distinguishing GRBs with typical photon energies below or above 0.8 MeV. 
We call them 'soft' and 'hard' bursts, respectively. We also carried out calculations for 
different choices of the critical c' range 2.5 <, c ~ ^ ^^^^ would correspond to 

different peak energies (~ 400 keV and 2 MeV, respectively; see Fig. 3 of T97). Our main 
conclusions do not depend on the choice of H^2,c range given above. Our preferred 

choice {HRl2^c — 3) is physically justified by selecting GRBs whose typical is ~ 4 times 
the average Ep detected by BATSE (~ 200 keV, Band et al. 1993), and outside the optimal 
energy range for BATSE detection (50-300 keV). The energy spectrum of GRB 910503 
Schaefer et al. 1996, T96a) clearly shows the spectral properties of GRBs corresponding to 
our critical value for H^2,c- '^^^ range for ideal BATSE detection is 1 <, H^2 ~ 2.5. 

Fig. 1 clearly shows that short GRBs have average Ep well above 1 MeV, and are 
efficiently detected by a combination of short risetimes (below 1 s) and spectral evolution 
in the BATSE's optimal detection range. We notice the relative small number of short/soft 
GRBs, contrary to the large number of long/soft bursts. This feature is not influenced by 
selection effects, and reflects the GRB hard-to-soft spectral evolution from simple GRBs 
to multiple-pulsed events. Selection effects may play a role in determining the trigger 
efficiency for long/hard events populating the upper right corner of the hardness/duration 
diagram. In our analysis, we consider only bursts with well defined spectral properties, 
including only events with relative uncertainty on H^2 ^^^s than 50%. We also note that 
only ~ 3% of the events in Fig. 1 show a deviation of their measured H^2 by more than 3a 
larger than the ideal maximum allowed SSM value {H^2 — 5, see T97). These events are 
most likely characterized by soft photon depletion/absorption in the 50-300 keV band that 
may influence the observed H^2 addition to instrumental effects. 
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3. Analysis of the 4B catalog 

Equipped witli a physical interpretation of the H^2 ~ ^90 diagram, we analyzed 
the brightness distribution of 4B GRBs according to their different hardness ratios and 
durations. We distinguish four classes: 

(S/H) GRBs with T90 < 2.5 sec, and H^2 > 3 (short/hard bursts); 
(S/S) GRBs with Tqq < 2.5 sec, and < 3 (short/soft bursts); 
(L/H) GRBs with T90 > 2.5 sec, and > 3 (long/hard bursts); 
(L/S) GRBs with T90 > 2.5 sec, and H^2 < 3 (long/soft bursts). 

Fig. 2 shows the results, where we plot the cumulative burst number with peak photon 
intensity above a given value Pq^ as a function of P^^ (64 ms trigger timescale). Both classes 
L/S and S/H show brightness distributions in good agreement with the —3/2 Euclidean 
slope. The deviations from the Euclidean distribution are minor, with a discrepancy below 
a factor of 2 between real and inferred values from the high-P64's at Pq4 = 1.3 ph. cm~^ s~^. 
The curvature near Pq4 1 ph. cm~^ s~^ is clearly due to a selection effect. Both L/S 
and S/H distributions might also be marginally consistent with a high-fiux slope of -3/2 
followed by an intermediate- flux slope near -1 for Pq4 8 ph. cm^^ s^^. On the contrary, 
the long/hard (L/H) bursts show a strong deviation from the Euclidean distribution at 
Pq4 ~ 20 ph. cm~^ s~^. The fourth class of short/soft events contains very few (66) and faint 
events and its brightness distribution is not relevant hereQ. These results are confirmed also 
for brightness distributions obtained for 256 and 1024 ms trigger timescales. 

All four duration/hardness GRB classes show isotropic sky distributions. Table 1 gives 
the values obtained for the dipole and quadrupole moments relative to the Galactic center 

^ If we summed together all the soft and hard events with T90 < 2.5 s we would obtain 
results very similar to those for the S/H class. 
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and plane. We also obtained the values of the V/V^ax parameter for GRBs with tabulated 
ratios Cp/ Cmin (with Cp and Cmin the peak and threshold countrates, respectively for the 
64 ms trigger timescale in the BATSE's energy range). Table 1 reports the average values 
of V/Vjnax for the three duration/hardness GRB classes. We note that a spectral bias 
in the determination of C^m is to be expected (see below), and that tabulated values of 
Cmin can be lower limits to the actual minimum countrates for detection of bursts with 
duration/spectral properties outside BATSE's optimal range. 

Our results indicate that only long/hard GRBs shows a strong deviation 
from the Euclidean brigthness distribution. These events dominate the bright 
end of the GRB population (one out of two GRBs are L/H bursts in the range 
6ph. cm^^s"^ < Pe4 <, lOOph. cm"^ s"^) and yet they show the strongest deviation from 
the Euclidean slope. 

4. Discussion 

We briefly discuss here some points relevant to the interpretation of our results, 
postponing a more thorough discussion to a later publication. Our results confirm 
previous analyses of the GRB brightness distribution, and contradict the predictions of 
simple cosmological models with no strong spectral and luminosity evolution. The total 
sample of GRBs detected by BATSE shows a clear deviation at low peak intensities from 
the Euclidean brightness distribution. This property, together with an average V/Vmax 
value near 0.37 (see Table 1) was interpreted as a signature of a globally inhomogeneous 
distribution in space (e.g., Paczynski 1991, Briggs 1996). However, we find that this 
conclusion can be applied only to the special subclass of long/hard bursts. The fact that 
only long/hard GRBs are strongly non-Euclidean being at the same time very intense is 
puzzling for extragalactic models. The soft and long GRBs that would naturally exhibit 
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an inhomogeneous brightness distribution in simple cosmological models show little, if any, 
deviations from being Euclidean. We emphasize that our conclusion is deduced by the 
combination of average spectral and time dilation properties of GRBs. Furthermore, the 
peak intensity distributions of S/H, L/S and L/H bursts are similar to each other, making 
unnatural that L/H bursts are in the foreground of more distant (and possibly softer) 
bursts. We are left with the conclusion that strong luminosity and spectral cosmological 
evolution is necessary to reconcile our results with extragalactic models of GRBs. L/H 
bursts placed in remote regions where a relative decrease of their spatial density makes 
them strongly non-Euclidean would have to be on the average very bright (by a factor 
of several) and with harder spectra {Ep larger by a factor of ~ 10) than the majority of 
closer bursts. Whether this can be accomplished by GRB formation mechanisms in star 
forming galaxies (SFGs) (as suggested by recent HST observations of GRB 970508, Pian et 
al. 1998) remains to be seen. The redshift dependence of SFG density (Madau et al. 1996) 
is strongly peaked near 2; ~ 1 and then decreases for larger redshifts. L/H bursts might 
originate at redshifts z ^ 1 under the condition that more distant GRB sources produce 
bursts more luminous and harder than those at smaller redshifts. The V/Vmax distribution 
of L/H bursts interpreted at face value with no selection bias, indicate that these bursts 
tend to reside in relatively close regions. The May 8th 1997 burst detected by BSAX (Piro 
et al. 1997) has an equivalent hardness H^2 — 3.6 ±0.9 and Tgo — Ss (Amati et al. 1997). If 
GRB 970508 belongs to the L/H subclass of GRBs, is in the strongly non- Euclidean part of 
the brightness distribution. Our results imply that, if the optical transient at cosmological 
distance is associated with the burst, GRB 970508 might be located at a redshift larger 
than z — 0.83 (indicated by optical spectroscopy, Metzger et al. 1997). 

Alternative interpretations of our results should also be considered. One possibility 
that reconciles our results with those of Pendleton et al. (1997) is that BATSE selection 
effects play a major role in shaping the brightness distribution of duration/hardness 
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subclasses of GRBs. The ability by BATSE to trigger on L/H events might crucially 
depends on their peak intensity (see also Lloyd & Petrosian 1997). As clearly indicated by 
the P97 results, even single pulses within complex GRBs share the same characteristics 
of our burst-averaged S/H, L/S and L/H subclasses. This fact demonstrates that short 
and long/soft bursts (pulses) can be detected with good efficiency given BATSE's trigger 
conditions. However, detection of relatively long and hard bursts (or L/H pulses usually 
occurring at the beginning of complex GRBs, see P97) require special conditions. It is 
conceivable that these trigger conditions are favorably met for quite bright GRBs above 
Pu ^ lOph. cm-^s"-*^. Detecting a burst (pulse) which is long/hard (or with substantial 
fiux above 300 keV) may be difficult for smaller values of Pq4. This might explain the sharp 
deviation from the Euclidean slope of the L/H burst brightness distribution. Since BATSE 
detects short/hard bursts with no apparent trigger inefficiency for low peak intensities (see 
Fig. 2), we conclude that long duration bursts (possibly with slow risetimes) with typical 
Ep^s outside BATSE's optimal detection range are the most affected by the trigger criteria. 
The results of Lloyd & Petrosian (1997) indicating a corrected Ep distribution with a large 
undetected population of GRBs with Ep 1 MeV agree with this interpretation. BATSE's 
trigger conditions [flux enhancement over background nominally by 5.5 standard deviations 
within a given timescale (64, 256 and 1024 ms) and for the energy range 50-300 keV] 
strongly favors detection of short bursts as well as long bursts with typical emission in the 
optimal energy range (our L/S bursts). Slow risetime-bursts with average Ep > 1 MeV 
and no strong hard-to-soft spectral evolution are disfavored for BATSE's detection. Both 
burst-averaged and individual-pulse properties of GRB brightness distributions share the 
same Euclidean (non-Euclidean) character of their soft (long/hard) population. This fact 
finds a natural interpretation in BATSE's selection effects. 

We also note that BATSE spectral selection effects can strongly influence the actual 
values of Cmin- Taken at face values, the average V/Vmax^^ obtained for our duration/spectral 
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GRB classes (Table 1) indicate inhomogeneous distributions by more than three standard 
deviations for all except the S/S bursts. However, BATSE's trigger conditions can couple 
the distance and spectral dependence of countrates in a non-trivial way, breaking the 1/r^ 
dependence of Cp deduced in previous apphcations of the V/V^ax method to GRBs (e.g., 
Schmidt, Higdon & Hueter, 1988). The nominal values of C^m tabulated in the 50-300 keV 
range can underestimate (by a factor of a few) the actual values of C^m necessary to detect 
bursts well outside BATSE's optimal spectral and risetime ranges. The long/hard and 
short/hard GRB populations are the most affected by this bias. 

Disregarding a Galactic origin for GRBs might have been premature. The combination 
of sky isotropy, inhomogeneity and the average value of V/Vmax < 1/2 apphed to the whole 
sample of GRBs was interpreted as being inconsistent with standard Galactic models (e.g., 
Paczynski 1991). However, our results clearly indicate that the strong deviation from a 
homogeneous distribution (if that is the correct interpretation of non-Euclidean GRBs) 
applies only to long/hard bursts, and certainly not to the majority of GRBs. A Galactic 
origin of a subclass or the majority of GRBs is not excluded by our results. 

The author thanks an anonymous referee for helpful comments. 
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Table 1: Properties of duration/ heir dness GRB subclasses 



GRB duration/hardness class 


< cos(^) > < sm^b) - 1/3 > < V/Vmax > 


long/soft (520) 
long/hard (338) 
short/hard (254) 
short/soft (66) 
total (1178) 


-0.002 ±0.043 -0.007 ±0.013 0.434 ± 0.019 (221) 
-0.020 ±0.054 +0.008 ±0.016 0.289 ± 0.022 (168) 
-0.062 ± 0.063 -0.019 ± 0.018 0.354 ± 0.023 (149) 
±0.003 ± 0.123 ±0.042 ± 0.036 0.458 ± 0.044 (43) 
-0.019 ±0.016 -0.007 ±0.013 0.373 ± 0.012 (581) 



The table gives the dipole, quadrupole moments with respect to the Galactic center 
and plane, and the average value of V/Vmax- We follow the definitions of Briggs et al. 
(1996). The angle 9 is the angle between the burst position and the Galactic center. The 
angle b is the Galactic latitude. 
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Fig. 1. — GRB energy fluence-averaged hardness ratio {H^2) duration (Tgo) distribution 
from the 4th BATSE catalog (4B). Only GRBs with fractional error on H^2 l^ss than 50% 
are plotted. Error bars are Icr. Long/hard (L/H) GRBs are marked by solid filled circles: 
they show the strongest deviation from the Euclidean brightness distribution (see Fig. 2). 
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Fig. 2. — Brightness distribution (cumulative number of bursts above a certain peak 
photon intensity level Pq4 as a function of Pq4 (64 ms trigger timescale) for different subclasses 
of GRBs. Solid histogram: total distribution; short- dashed histogram: long/soft bursts; long- 
dashed histogram: long/hard bursts; dotted-dashed histogram: short/hard bursts. 



